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m Abstract 

The r-mode instability in rotating compact stars is used to constrain the phase of matter at high density. 
The color-fiavor-locked phase with kaon condensation (CFL-KO) and without (CFL) is considered in the 
temperature range 10 8 K < T < 10 n K. While the bulk viscosity in either phase is only effective at damp- 
ing the r-mode at temperatures T > 10 n K, the shear viscosity in the CFL-KO phase is the only effective 
damping agent all the way down to temperatures T > 10 8 K characteristic of cooling neutron stars. How- 
ever, it cannot keep the star from becoming unstable to gravitational wave emission for rotation frequencies 
V « 56 — 1 1Hz at T « 10 8 — 10 9 K. Stars composed almost entirely of CFL or CFL-KO matter are ruled out 
by observation of rapidly rotating neutron stars, indicating that dissipation at the quark-hadron interface or 
nuclear crust interface must play a key role in damping the instability. 
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I. INTRODUCTION 

The r-mode instability of rotating compact stars is of wide astrophysical interest for several 
reasons [1 J. The majority of neutron stars have rotation period P ~ 0.5seconds, much slower than 
theoretical limits: P ~ 0.001s (v = \/P ~ 1kHz) [2]. Even the fastest spinning neutron star PSR 
J1748-2446ad with v = 716Hz is not near the theoretical stability limit, with the majority of 
"spun-up" neutron stars in low mass X-ray binaries (LMXBs) catalogued between 300-640 Hz. 
The r-mode provides a mechanism for spinning down young neutron stars [4J and limiting spin-up 
frequencies in millisecond pulsars (see eg. Q for a review) through angular momentum loss in 
gravitational wave emissions. These gravitational waves may be detected given the expected im- 
proved sensitivities of ground-based interferometers eg. VIRGO, advanced LIGO [6|. In addition, 
r-modes are probes of the equation of state of cold and dense matter through viscous damping 
effects and advance our understanding of nuclear and quark interactions by providing an astro- 
physical and phenomenological link to Quantum Chromodynamics (QCD) over a certain range 
of temperatures (10 7 K < T < 10 U K) and at super-saturation densities [rig ~ (2-5)n sa t]. Only a 
few studies of r-modes linking them to phenomenology of quark matter in neutron stars have been 
carried out thus far [714111. In this work we use r-mode damping to constrain the presence of a 
realistic phase of superconducting quark matter lfT2l[P3l in the core of neutron stars: viz., one that 
includes the effects of a sizable strange quark mass m s at intermediate baryon density. 

A key parameter in understanding the r-mode instability in rotating compact stars is the critical 
frequency Q. c for the onset of the instability (Q. = 2n/P rad/s). Gravitational radiation tends to 
make the r-mode grow on a characteristic time scale Xqr- Internal friction such as viscosity or 
mutual friction tends to damp r-mode growth on a characteristic time scale Xp. The competition 
between the two determines Q, c at any temperature. At Q. > Q. c , the r-mode develops an instability 
(1/ttotai = 1/xqr + I/^f < 0). Previous works [|9]-[TTI have studied the r-mode damping in normal 
(ungapped) quark matter as well as some superconducting (gapped) phases such as the color- 
flavor-locked phase (CFL), and found that bulk and shear viscosity in the fully gapped phase 
is large enough to render rapidly spinning quark stars (£1 > 0.1 ^Kepler) stable. However, this 
conclusion applies only in a restricted temperature regime T > 10 10 K, since the shearing mean free 
path due to phonons otherwise becomes larger than the star's radius (~ 10km). Furthermore, for 
temperatures T < 10 8 K, mutual friction in the CFL phase has been shown to be too weak to damp 
r-modes |fTT|. arguing for rapid spin-down of CFL stars to less than a few Hz, effectively ruling 
our cold CFL matter in any rapidly rotating neutron stars. The r-mode damping in the temperature 
regime 10 8 K < T < 10 10 K for color- superconducting phases of quark matter is therefore an open 
question, and is the focus of this work. Here we study r-mode damping and the critical frequency 
curve for a neutron star made up mostly or entirely of the kaon condensed CFL (CFL-K0) phase. 

The CFL phase with symmetric pairing of up-, down- and strange-quarks is severely stressed at 
realistic chemical potential /j q ~ 300 MeV because the strange quark mass m s ~ 100 MeV is non- 
negligible and costs extra energy m^/(2fj q ) compared to the much lighter up- and down-quarks. 
In Ref. Ifl4l . it was shown that it is energetically favorable for the CFL vacuum to be in a rotated 
chiral state with a non-zero kaon condensate when the cost m^/(2/u q ) exceeds the kaon mass m^, 
the lightest meson in high-density QCD [fT51 . Realistic neutron star densities could support the 
CFL-K0 phase. We find that in the CFL-K0 phase, while the bulk viscosity has little role to play 
in the temperature region relevant to r-mode damping in cooling neutron stars, the shear viscosity 
from kaons in the CFL-K0 phase is important even below T < 10 9 K, unlike the CFL phase. In 
essence, this fact coupled with the large mean free path of the phonons controls the main features 
of the critical frequency curve (Fig. [2]). 



II. THE r-MODE 



r-mode is a pulsation mode intimately linked to the rotational properties of the star, and the 
restoring force here is the Coriolis force. The frequency of r-mode oscillations depends on the 
star's rotation frequency as well as the average density of the star. Including second-order rota- 
tional effects, the mode frequency co r in the inertial frame is |[Tl[T6l 
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where oo rot is the mode frequency in the co-rotating frame and m is the azimuthal quantum number 
of the mode, taken to be 2 for the first r-mode to become unstable (higher m modes become 
unstable at higher frequency). The number K2 includes the sensitivity to the density profile (see 
Refs. [|9l[I3), po is the average density of the unperturbed star and G is Newton's constant. 

The r-mode couples to gravitational waves through the current multipole of the perturbation ||T). 
All modes with azimuthal quantum number m > 2 suffer the so-called CFS instability |[T8l [T9ll 
when £1 > (H rot Jm so that the r-mode energy grows at the expense of the star's rotational kinetic 
energy as angular momentum is lost to gravitational wave emission. In general, viscous forces 
counter this energy growth. The timescale x associated with growth or dissipation is given by 
\/%i = —jg {jzn ■■> where the subscript i corresponds to various dissipative forces. Details about 
computing E and (jjr) ■ are given in Refs. []9l HVj . Explicitly, the gravitational radiation timescale 
is KB 
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This timescale is negative, indicating exponential mode growth (the instability). 

Viscosity (bulk, shear and otherwise) is a source of damping for the r-mode. For a given phase 
with shear viscosity T], the shear viscosity timescale is Q]| 
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We find that the bulk viscosity damping timescale plays a minor role in our calculation, see Table|I] 
and Figs. [T}[2} 

Mutual friction is another source of dissipation for the r-mode energy. It occurs due to the 
scattering of phonons from quantized vortices in the superfluid component. For the CFL phase, 

this has been calculated to be [fTTI : ^— ~ 36.2 ( — J Q. in the limit of negligible vortex mass. The 

critical rotation frequency Q. c of the star can be determined by the criterion that at this frequency, 
the fraction of energy dissipated per unit time exactly cancels the fraction of energy fed into the 
r-mode by gravitational wave emission: 
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Stable rotation frequencies at any temperature are limited by the smaller of the critical fre- 
quency as computed from Eq. ^ or the Kepler limit Q,/Q,k=1. As depicted in Fig. [2j the region 
above the temperature-dependent Q. c curve is unstable to r-mode oscillations and the star upon en- 
tering this region will be spun down rapidly to Q. < Q. c by emitting gravitational waves and losing 
its angular momentum in the process. 



III. CFL-KO: EOS AND VISCOSITY 

For simplicity, we adopt a homogeneous composition assuming that the entire star is made of 
CFL-KO matter. This allows a direct comparison with r-mode calculations in hadronic matter with 
and without a crust [|20l . We will see later that a crust is required for agreement with observations, 
but the more elementary case of a pure CFL-KO star should be tested first. Just as for the CFL 
phase, the neutral kaons are expected to dominate over the phonons in determining the bulk viscos- 
ity. However, the condensation changes the temperature dependence of the viscosity, which is the 
most important outcome for the r-mode damping. The pressure of the CFL-KO phase, including 
terms up to 0{m A s ) is given by lfl4ll2TI 
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The pressure from the Nambu-Goldstone bosons introduces a negligible temperature dependence. 
We use the perturbative expression for f n |fT31 . This EoS differs from the CFL EoS only at 0(m A s ) . 
In the above expression, the common Fermi momenta kp also determines the CFL pressure Il22ll 
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where A is the superconducting quark gap and B is the MIT Bag constant. We choose A = 100 MeV 
and B = 80 MeV/fm 3 such that the energy per baryon E/A < 930 MeV at zero pressure. In the 



limit of vanishing up and down quark mass, we have up to 0(m A s ): kp = 

/j-m 2 /(6/u)+m*/(72v 3 ). 

Evaluating the expression for Pcfl-ko pressure using Eq.Q, we obtain 
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where f 2 = f^/kp pa (21 — 81n2)/(367t 2 ) . The energy density then follows from standard thermo- 
dynamic relation £ = — Pcfl-ko + n/u with number density n = 3 /U /'cfl-ko- Eliminating /j between 
pressure and energy density, we get the EoS in the useful form 
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The dominant contribution to the bulk viscosity in the fully gapped CFL or CFL-KO phase is 
determined by weak reactions involving the lightest modes Il23ll24l . In the CFL-KO phase, this is 
the neutral kaons, a fraction of which are in the condensate. From Alford et. al ll24l . bulk viscosity 
is 
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where the kaon number density n, susceptibility % and weak equilibration rate for kaons in the 
condensate are computed in Ref. Il24l using the 2PI formalism in order to assess thermal effects 
consistently. Note that dm = m# —/jk<0 for the condensedphase but dm > for the pure CFL 
phase. For the typical parameter values flU |23] used in Fig. [TJ T c ~ 3 x 10 11 — 4 x 10 11 K. There 
is a negligible contribution to bulk viscosity due to phonon scattering above T > 5 x 10 8 K fl25). 

Kaon condensation breaks flavor symmetry, and the associated Nambu-Goldstone mode con- 
tributes to shear viscosity at low temperature. Although the kaonic contribution is smaller than 
that from the phonons, it becomes important when the phonon mean free path exceeds the stellar 
radius at T < 10 10 K - see Fig[T] The leading order contribution was calculated in terms of an 
unknown dimensionless coefficient C [|26l . Analytic expressions for the shear viscosity could be 
derived for very small or large value of the coupling C: 

11 Ma4 
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For natural values C ~ 1, shear viscosity r\ has to be calculated numerically which we have done 
and cross-checked against the results of Ref. Il26l . Here also r\ shows a T 5 dependence similar 
to the contribution from massless phonon scattering. Phonon contribution to shear viscosity is 

riPhonon = 3.745 x 10 6 (^) r- 5 g/(cm s) , (11) 

where Tg is the temperature in units of 10 9 K. The effects of mutual friction have been shown to be 
too small to effectively damp the r-mode instability at T < 10 8 K [[TT]|. In this work, we include 
the mutual friction damping. 

In Fig. [T] below, we show the bulk and shear viscosity from the CFL-KO phase in comparison 
to the pure CFL phase. 



IV. DAMPING EFFECTS 

Next, we present results for the damping times from the various mechanisms discussed pre- 
viously and obtain the critical frequency plot. The damping times are listed in Table [TJ The bulk 

TABLE I. Damping times at T = 10 9 K for the CFL-KO phase, the CFL phase with 0{m A s ) corrections 
(labeled CFL4 ) and without. For the CFL-KO phase, the shear damping timescale due to kaons with 
parameter value C = 1 (see text) is shown. The central energy density p c (in units of nuclear saturation 
density po = 2.54 x 10 14 g/cm 3 ) is chosen to yield a 1.4 M star. The stellar rotation frequency is taken to 
be half the Kepler frequency Q.k = f \/2nGp. 

EoS p c /po j?(km) Q. K (rad/s) x^ (s) x^ (s) x G w (s) 

CFL-KO 2.526 10.748 6696 9.6 xlO 22 4.9 xlO 13 -7.6 xlO 2 

CFL4 2.526 10.709 6697 1.7 xlO 12 - -7.8xl0 2 

CFL 2.526 10.729 6696 1.7 xlO 12 - -7.7 xlO 2 
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FIG. 1. Top panel: viscosities for the CFL and CFL-KO phase. Bottom panel: phonon and kaon shear 
viscosity mean free path X. The typical parameter values |9||22| are fi q = 310 MeV, /Uk = 17.92 MeV and 
CO = 271 ms _1 , 8m = 1 MeV for CFL, 8m = — 1 MeV for CFL-KO. In the bottom panel, the crossings with 
the solid (black) horizontal lines at 1km and 10km indicate when X is on the order of the stellar size, making 
the shear viscosity irrelevant. 



viscosity of the CFL-KO phase is smaller than that of the CFL or ungapped phase for temperatures 
T < 10 K. Consequently, bulk viscosity damping timescales are much longer than the typical 
timescale for gravitational wave growth, rendering bulk viscosity irrelevant to r-mode damping in 
cooling neutron stars. The damping timescale from shear viscosity in the CFL-KO phase becomes 
small enough to effectively damp the r-mode as temperatures fall to T ~ 5 x 10 10 K. However, 
this is still dominated by the phonons, since the kaon shear viscosity is much smaller. The phonon 
mean free path increases with decreasing T, therefore, phonon shearing collisions occur less often. 
The mean free path X§ is estimated from the kinetic theory relation X§ ~ r\/(n§p§) using phonon 
density n§ and thermal momentum p§, see Refs. D9l|26l|271. In Fig.fTl we show the mean free path 
of the phonons and kaons. A larger value of the parameter C corresponds to a larger cross-section 
for shearing collisions and therefore yields a smaller mean free path at the same temperature. It 
can be seen from Fig. [Tjthat kaon shear viscosity is important in the temperature range 10 8 K 
< T < 10 10 K, where phonon shear viscosity is no longer important. 

Finally, we compute the critical frequency curve for the CFL-KO phase according to Eq.Q. 
We notice two "dips" in the critical frequency curve: the first one at T ~ 8 x 10 9 K, when the 
phonon mean free path X§ becomes larger than the stellar radius, and the second at T ~ 2.5 x 10 8 
when the kaon mean free path becomes larger. This second dip is sensitive to the value of C, 
appearing earlier for smaller C. In order to obtain a smooth behaviour for the critical frequency 
curve, we use for convenience a sharply peaked cutoff function f(kty) = [1 — tanh^/km — 5)]/2 
centered around a 5 km stellar radial distance. The maximum stable frequency accommodated by 
the CFL-KO phase falls well below the observed bound on rotation rates in LMXBs. 
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FIG. 2. The critical frequency v as a function of temperature for CFL-KO quark stars. The curves end 
abruptly when the shear mean free path exceeds the stellar radius. We use the same parameters as in Fig.[T] 
The box represents typical temperatures [28] (2 x 10 7 -3 x 10 8 K) and rotation frequencies (300-700 Hz) of 
the majority of observed neutron stars in LMXBs. 



CONCLUSIONS 



r-mode oscillations of compact stars in the kaon-condensed CFL phase are considered. The 
mode frequency is almost exactly the same as that for a CFL star [9] since the equation of state for 
the CFL-KO phase only differs at O(m^). The mode growth and viscous damping timescales in the 
temperature range 10 8 K < T < 10 n K is studied based on the dominant bulk and shear viscosity 
contributions. Compared to the pure CFL phase, the bulk viscosity in CFL-KO is smaller and does 
not play a significant role. At temperatures T > 10 10 K shear viscosity associated with phonon scat- 
tering is important and determines the temperature dependence of the critical frequency curve. At 
lower temperatures, the phonon mean free path exceeds the star radius, so that the new mechanism 
for shear viscosity associated with kaon condensation 11261 provides the dominant dissipation. In 
comparison to pure CFL stars, CFL-KO stars are more stable against the r-mode instability at the 
lower end of the temperature range considered. However, given the observed LMXB spin rates, the 
critical frequency curve predicts that even pure CFL-KO stars are unlikely to exist. Just as in the 
neutron matter case, viscous damping just beneath the core-crust interface (Ekman layer) appears 
to be necessary to provide a large enough stable rotation rate consistent with LMXB data [20] . 

The main result of the current work is that given the known dissipative mechanism in CFL or 
CFL-KO phase, a pure quark star with these phases is ruled out by observed LMXB spin rates. Mu- 
tual friction associated with kaon-vortex scattering could provide an additional dissipative mech- 
anism, but the critical frequency curves calculated here provide a lower bound. Future related 
work should investigate hybrid stars with a CFL-KO core. The role of the crust-core interface as 
well as the quark-hadronic matter layer should be key ingredients in obtaining constraints on the 
persistence of gravitational waves from stars containing superfiuid quark matter in their interior. 
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